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Abstract: We consider a brane generated by a scalar field domain wall configuration in 
4+1 dimensions, interpolating, in most cases, between two vacua of the field. We study the 
cosmology of such a system in the cases where the effective four-dimensional brane metric 
is de Sitter or anti de Sitter, including a discussion of the bulk coordinate singularities 
present in the de-Sitter case. We demonstrate that a scalar field kink configuration can 
support a brane with dS4 cosmology, despite the presence of coordinate singularities in the 
metric. We examine the trapping of fermion fields on the domain wall for nontrivial brane 
cosmology. 
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1. Introduction 

Models incorporating more than 3 + 1 dimensions have been of great interest to the high- 
energy physics community from the early 1980's onward, and have been investigated in 
considerable depth over the last decade (see, for instance, Refs. ||l], ^, Q). It was initially 
thought that any such extra dimensions would need to be compactified to a small radius in 
order for their effect on four-dimensional gravity to be consistent with current experiments. 
However, the work of Randall and Sundrum showed that effective four-dimensional 

gravity could be recovered even in the case of non-compact extra dimensions. 

The original Randall-Sundrum proposals, and much of the subsequent work, dealt 
with Standard Model fields confined a priori to a 3+1-dimensional subspace of the higher- 
dimensional space, with gravity confined to this "brane" by a warped five-dimensional 
metric (reviews are given in |5|, P). In these models the brane is simply the boundary 
between two domains in the bulk spacetime, with the energy-momentum localised on the 
brane corresponding to a discontinuity in the derivatives of the metric across the brane. 
The ensuing cosmological evolution of the brane-bulk system has been studied extensively 
in the literature (|^, or ^ for discussions of cosmology induced by moving branes). 
Many modifications and extensions to the initial RS scheme have also been proposed: see, 
for example |, |, |l|, 0, |l3|. 



It has been suggested that such a 3+1-dimensional brane could also be realised as a 
smooth domain wall that dynamically confines the fields of the Standard Model (|l], 15 
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P^ , |r^ , |T^, |l^). General aspects of smooth analogues of the Randall-Sundrum model have 



been studied in ||2^, 21 1. In particular, given a scalar field potential with two neighbouring 
minima, a domain wall soliton may form where the scalar field interpolates between two 
regions of spacetime corresponding asymptotically to vacua of the field. This approach is 
appealing because it provides a dynamical origin for a brane model with one infinite extra 
dimension. 

A number of authors have investigated solutions to the Einstein equations for a five- 



dimensional warped (nonfactorisable) metric coupled to a scalar field 15, 16, |l^, |l9| 



22, 23, 24, 25, 26, see also ||2^, for non- minimally coupled scalar-gravity models). 
The trapping of fermion fields on a subspace in models with extra dimensions has also 
been explored ( |l], |l|, |l^ |l|, ||, ||, ||, ||]). However, most (but not ah) of these 



investigations have only considered the case where the effective 4D metric is Minkowskian, 
so they do not explore the cosmology of these models. The purpose of this paper is to 
study aspects of the cosmology of domain-wall-style branes. 

In Sec. |2| we write down the Einstein equations for a warped 5D metric based on an 
arbitrary 4D metric, with the energy-momentum tensor given by the bulk scalar field. We 
then discuss the properties of solutions to the equations of motion, in the cases where the 
effective 4D metric is Minkowski, dS4 and AdS4. 

After reviewing some examples of scalar field domain walls with 4D Minkowski slices 
(Sec. we present criteria for a smooth 5D metric containing no curvature singularities 
with dS4 slices in Sec. ^ In Sec. ^ we write down one such metric analytically, and plot 
the corresponding scalar field and potential for a particular parameter set. In Sec. ^ we 
present new analytic solutions for the 5D warped metric with AdS4 and dS4 slices on the 
brane, and the associated background scalar field, motivated by thin-brane solutions that 
have been employed to demonstrate localisation of gravity ( ||36[| ). 

Finally, in Sec. ^ we consider the trapping of fermions in a spacetime described by a 
5D warped metric (with dS4 or AdS4 brane cosmology) and supported by a kink-like scalar 
field domain wall. In Sec. |8| we examine the behaviour of the chiral fermion zero modes in 
these various solution systems. 



2. The equations of motion for a warped metric coupled to a scalar 

Suppose gj^J 
simple form. 



Suppose gl^J is some general 4D metric. Let us consider a 5D warped metric guN of the 



mN=[' 7- J). (2.1) 

(In general, we shall use capitalised Roman letters to indicate five-dimensional indices, 
and Greek letters for four-dimensional indices.) We use the sign conventions of |37]: in 
particular, the 5D metric has signature ( — h + + +). 

This is a generalisation of the usual Randall-Sundrum warped metric ansatz ( @ ) to 
the case of a non-trivial four-dimensional metric. Such cosmological solutions have been 
studied in the context of infinitely thin branes by a number of authors (||6|, 0, ^, |l^). We 
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shall refer to e"'^*-^^ as the "warp factor" associated with the warped metric, and cr{y) as 
the "warp factor exponent". 

This simple ansatz clearly does not encompass all possible 5D metrics that yield effec- 
tive 4D gravity. In particular, if the bulk component T44 of the energy-momentum tensor 
depends only on the bulk coordinate, then the Einstein equations require that the effective 
4D metric generated by this ansatz must be of constant curvature, allowing dS4 and AdS4 
brane cosmology, but prohibiting more general FRW solutions. 

It is well known that in the case of an infinitely thin brane, the presence of matter on 
the brane causes mixing of the bulk and brane coordinates in the metric, in the coordinate 
system where the brane is a hypersurface of constant y (||8|, ^). This effect is frequently 
termed "brane bending". A simple concrete example is shown in Q, where an FRW 
ansatz is used for the 5D metric, the brane lies at y = 0, and the metric components are 
separable into bulk and brane coordinates only so long as the brane tension is the sole 



contributor to the energy-momentum tensor. However, the ansatz of Eq. |2.1| is useful 
for studying systems where the brane tension dominates the matter contribution to the 
energy-momentum tensor, as in some cases the ensuing Einstein equations can be solved 
exactly. 

Denoting the Ricci scalar computed from the 5D metric gM N by R^^^ , and the 4D Ricci 
scalar computed from g^J by we find that the 5D Ricci scalar and Einstein tensor 
can be written in terms of their 4D counterparts and the warp factor as per 

(G'^^U3gi!Je-My)(a"iy)-2{a'{y)f) \ 

""""-I, -leM.)^(4)_6(,^(y))^j' 

i?(5) = 62^^(2^) + (20 - 8a" (y)) . (2.3) 

The Einstein equations in five dimensions are. 



Gmn = -SttGTmn, (2.4) 
and in particular, writing 8ttG = k^, we see that 

^^r44=QeM^)iZ(^)+6(a'(y))^), (2.5) 

implying that a time-varying 4D curvature scalar is associated with a time-dependent r44 
(similarly, variation of the curvature scalar with the other brane coordinates requires r44 
to depend on those coordinates). Thus if the brane tension is independent of the brane 
coordinates {t ^), and the trapped matter (if present) does not contribute to r44, 

only constant-curvature solutions are possible with this metric ansatz. 

Since we expect these four-dimensional metrics to be valid only in the case where the 
brane tension dominates over the matter contribution to the energy-momentum tensor, let 
us compute the 5D energy-momentum tensor derived from a scalar field ^ and associated 
potential V{^). 
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In the absence of matter sources, the system is described by the usual Einstein-Hilbert 
action with additional scalar field terms ([^, |3^), 



S 



d^x\/ —g{x) 



2^2 



g{x) = Det {qmn) ■ 
The resulting energy-momentum tensor is easily computed as, 

■.MN 



(2.6) 

(2.7) 
(2.8) 



I =2 \- g 1,$ 

09MN 



V{'^) + ^g'^'^dR<^ds<^ 



(2.9) 



Suppose the scalar field $ is a function of the bulk coordinate y only. Then using the 
separable metric ansatz of Eq. 2.1, the non-zero elements of the energy-momentum tensor 
become. 



The Einstein equations, Gmn = —i^ Tmn, and Eqs. 2.2-2.3, then yield: 



^(4) 



-2<7{y) (4) 



(v{<^) + \ {^'{y)f) - 3 {a"{y) - 2 {a'{y)) 



^(4) ^ 2e-2-to) 



-6 {a'{y)Y + ^'(-{^'{y)Y-Vm 



(2.10) 
(2.11) 

(2.12) 
(2.13) 



19'^} (^ R^^^ = Then we obtain the equations of motion, 



For these equations to be consistent, there must be some constant 7 such that G 



(4) 



76 



My) 



^~\V{^) + -{^\y)y 1 -3 



a"{y)-2{a'{y)f), 



(2.14) 

= 3 {<y'{y)f - y Q {^\y)f - V^^ . (2.15) 

We can eliminate the potential from these equations and obtain the result, 

3a" (y) - {^'{y)f = ^e^^'^yy (2.16) 
The potential V can easily be obtained as a function of y from either of Eqs. 2.14- |2.1^ , 

^(y) = \ {^\y)f - ^ [o\y^^ + 5^'^^'^- (2.17) 

Zi r\j rv 

In the case of AdS4 space, 7 is negative, while for dS4 space the converse holds true. 
Minkowski space corresponds to the case 7 = 0. 
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For 7 7^ 0, we can define the dimensionless quantities, 

V = 

Hv) = 




In terms of these dimensionless quantities, Eqs. p.l6| and |2.17] become. 



g2a(r,) ^ M 

7 
1 



V{v) 



$'(7?) -2(a'(r/))' + 2 



l7l ■ 



(2.18) 
(2.19) 

(2.20) 

(2.21) 
(2.22) 



Note that, independent of the value of 7, differentiating Eq. 2.15 and then substitut- 
ing Equation p. 16 yields the Klein-Gordon equation for the scalar field (this relationship 
between the Einstein equations and the Klein-Gordon equation is discussed in pT 



dV 
d¥ 



^"{y)-Aa\y)<^>\y), 



or using the dimensionless quantities defined above, 

^ = |,"(^)-4a'(r?)<l'(,?). 



(2.23) 



(2.24) 



Later we will wish to discuss configurations where the metric elements go to zero at 
some point in the bulk. In this case, a is an unsuitable parameterisation as it diverges 
to 00 at zeroes of the warp factor. If in Eq. 2A we replace e"^"^ with w^, where u; is a 
real function of the bulk coordinate, then the equations of motion (Eqs. p. 16 -2. 17) can be 
computed exactly as previously, and become, 

3 (a;'(y))' - 3u;{yW'iy) - {^'{y)fu:{yf = 7, (2.25) 

2 



V{y) 



6 (J{y) 



+ 



27 



(2.26) 



\ u){y) J K?uj{y)'^ 

(Note: It immediately follows from these equations of motion that if the first and second 
derivatives of lo with respect to the bulk coordinate remain finite, then lo may possess 
zeroes only if 7 > 0, i.e. where the 4D metric is dS4 or Minkowski. We will later show that 
zeroes in lo inevitably arise for the case of a smooth oo with dS4 brane cosmology.) 

In terms of the dimensionless quantities defined in Eqs. 2.18 - 2.2C , the equations of 
motion for uo are, 

2 



171 
7 

V{ri) 



[UJ 



$'(7?) 



2 



uj{ri) ) 



+ 2 



7 



1 



I7I ^iv) 



2 ■ 



(2.27) 
(2.28) 
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We will preferentially use the a notation when not dealing directly with cases where 
cr — > oo, for ease of comparison with the literature. Coordinate singularities in the metric, 
where a diverges (and uj goes to zero), correspond to bulk horizons, and most previous 
braneworld studies have been concerned only with the behaviour of the brane-bulk system 
in the region bounded by these horizons. 

3. Scalar field kink configurations supporting a Minkowski brane 



In the case where 7 = and the brane metric corresponds to Minkowski space, Eq. p. 16 



can be integrated to give a warp factor corresponding to any chosen scalar field profile. The 



potential generating this scalar field can then be analytically determined from Eq. p. 17 
with 7 = 0. Many solutions to the coupled gravity-scalar system for a Minkowski brane 
have been discussed in the literature (see for example |l^, |l^, |2^, |3l| ) . 

Two simple examples are the kink configurations = (\/3^/k) tanh(ry), ^2(11) = 

(y/SA/ k) avctansmh{ry), where A is a dimensionless constant and r has dimensions of 
inverse length. Eq. |2.16| then yields the corresponding warp factor exponents 0"i,(T2: 

cri(y) = (-lncosh(ry) sech^ (ry) + Bry + C] , (3.1) 

V3 6 J 

(y) = A^ (In cosh(ry) + Bry + C) , (3.2) 

where B and C are dimensionless constants of integration. If the warp factor exponent is 
an even function of y (a usual symmetry of domain-wall solutions and the RS model with 
a single brane 0]), then it follows that the integration constant B must be zero. Then in 
both cases, the warp factors e~'^ converge to zero far from the brane, like the warp factor 
of [0] for the case of a fine-tuned infinitely thin brane. The corresponding potentials Vi, V2 



can easily be obtained from Eq. 2.17 : Vi is a sextic polynomial in <I> ( |18|| ), while 1^($) is 
3^2 r2 



2 «n('-'' + ''™'(i7i*)+'-''''™(i7i*) + 

+ {aA^B^ - 1)] . (3.3) 



If B vanishes (required if the warp factor is even) then the potential is also an even function 
of ^>. 

This method is useful for deriving analytically tractable solutions, but is conceptually 
somewhat unsatisfactory because the system is fundamentally determined by the poten- 
tial, not the profile of the scalar field. "Superpotential" methods have been used to reduce 
the equations of motion to a system of first-order equations by introducing an auxiliary 
potential (|l^, ^ 0); allowing the scalar field profile and metric corresponding to a wide 
range of potentials to be derived. These calculations can often be performed analytically, 
although writing the potential as a function of the scalar field rather than the bulk coordi- 
nate requires an inversion of the scalar field profile which is often analytically intractable. 

If 7 is non-zero, i.e. the 4D metric has non-zero constant curvature, then the problem 
becomes significantly more complex and this superpotential approach no longer yields easily 



solvable equations of motion. This problem has been investigated in |1£, 23, p^]) 
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A modified superpotential method can be employed to generate scalar field configura- 
tions and metric warp factors corresponding to a given auxiliary potential ([^]). However, 
it is applicable only to a limited class of potentials. The solutions found by this method 
with AdS4 brane cosmology exhibit kink-like scalar field configurations, but the dS4 solu- 
tions contain naked curvature singularities in the metric, associated with infinities in the 
scalar field as the potential becomes unbounded below, at a finite distance from the brane. 

4. Properties of solutions to the Einstein equations with dS4 cosmology 

Putting the superpotential approach aside, let us consider the general properties of solutions 
to the equations of motion (Eqs. |2. 21 -2. 22) when 7 > 0. We will prove that in this case, 



corresponding to dS4 brane cosmology, there are no smooth, even solutions for the exponent 
a{r]): if the warp factor e~'^^^^ is to be smooth and even, it must possess zeroes. 

We begin by noting that for 7 > 0, a"{r]) = (<l'(?7))^ -|- e^'^^''^ > OVr/. Thus cr'(r/) is 
monotonically increasing Vr/. Suppose that cr{r]) is an even continuous function of r] with 
continuous and well-defined first and second derivatives: then cr'{0) = and consequently 
cr'{ri) is negative for r/ < and positive for r/ > 0. Thus if cr^ri) is a smooth function of 
r], then it must be strictly decreasing for r] < and strictly increasing for r] > 0, and it 
follows that cr(0) is a minimum of the metric exponent a: cr{r]) > cr(0) Vt/ 7^ 0. 



Substituting back into Eq. |2.21| , it follows that 

cj"(77) >e2-('') >e2'^W, (4.1) 



and integrating this inequality implies that cr(r/), |cr'(r/)| — > cxd as 77 — > ±00. 
Let g{rj) = a'{'q)e~"^'^\ Then we can rewrite Eq. |2.21 in the form. 



By hypothesis, g{r]) is a smooth odd function of rj, with g(ri) > for r/ > 0, and ^'(0) = 

C7"(0)e-'^(0) y Q 

Now wherever g'{'r]) < 0, we must have g{r])'^ > 1, which for rj > implies g{r]) > 1. 
Taking the contrapositive, in the range rj > 0, whenever g{rj) < 1, g'{i]) > 0. Thus if 
g{i]) < 1, Vr/ > 0, then g is bounded and monotonically increasing and converges to some 
(strictly) positive real number S < 1. Furthermore, since g{rj) is monotonically increasing 
on any interval for which g^rj) < 1, and g is continuous, it follows that if 3r]Q s.t. g{rjQ) > 1, 
then T] > TjQ ^ g{r]) > 1. 

In both of these cases, 3 some real and strictly positive number e, and some rjQ > 0, 
such that rj > r]o ^ g{rj) > e. (In the first case, we can simply choose e = 6/2.) 

But then choosing some r/i > r/o, and integrating g{rf) over this interval, we find. 



VO -J VO 

(4.3) 

But this is a contradiction, as by choosing ryi arbitrarily large we can obtain an arbitrarily 
large lower bound on , which must be finite. 
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Consequently, there is no even, smooth warp factor exponent defined on the 

whole real line that yields a four-dimensional deSitter cosmology. It is possible that the 
warp factor io{y) is smooth, but in this case it must contain zeroes (as smoothness of 
CO = zbe"*^ implies smoothness of a except where to changes sign). 

It is well known that warped metrics with dS4 slices often contain horizons or naked 



singularities at a finite distance from the brane ([0, 19, 25, 26, ^3]). The properties of 



horizons in the case of a thin brane with no scalar field have been discussed in [^, 13 1 



and curvature singularities in the presence of a scalar field were considered in 25 1. We 
have demonstrated the strong result that such features are inevitable in the case of a brane 
supported by a scalar field. 

In the simplest analytic solutions for this system, the scalar field diverges to infinity 



at these points, and naked curvature singularities occur ([p3|, 24, ^]). We will show that 
it is possible to avoid a diverging curvature at these points, and to obtain a kink-like 
configuration for the scalar field supporting the dS4 brane, but at present we have only 
obtained a partially analytic example of such a solution. 

The physical interpretation of singularities in the metric, for a model where gravity 
is coupled to a scalar field, has been discussed in the literature ([^, 26|). In particular. 



Gremm suggests that five-dimensional spacetimes consisting of domain walls interpolating 
between spaces with naked curvature singularities may be interpreted as four-dimensional 



gravity coupled to a non-conformal field theory |25]. Gremm also claims to present a 
warped metric with dS4 slices where the curvature remains finite at the metric zeroes, but 
the scalar field and potential diverge, but we have not been able to verify this solution. 
Davidson and Mannheim give a similar solution in the case of an infinitesimally thin brane, 
with a divergent scalar field and potential but bounded curvature at the metric zeroes |Q , 
and suggest that such a solution might provide a mechanism for dynamical compactification 
of the extra dimension. 

To distinguish between curvature and coordinate singularities, we examine the five- 
dimensional curvature scalar, given by Eq. |2j. Also, -k^T^ = Gjfj = R^j - {l/2)gl^R = 
R{1 — n/2), where n is the number of dimensions. So the five-dimensional curvature scalar 
can be expressed in terms of the trace of the energy-momentum tensor, 

R^'^ = ^-^T^l- (4.4) 

In the case where the energy-momentum tensor is generated entirely by the scalar field and 
potential, 

Tji! = -4 (vi<^) + i (l>'(y))') + I {^'iy)f - V{<^) 

= -5m-^(^'(2/))'- (4.5) 

Thus, the curvature scalar can only become singular if the potential, the scalar field deriva- 
tive, or both, diverge to infinity. If the potential and scalar field remain bounded at a zero 
in the metric, this is sufficient (albeit not necessary, as diverging terms on the RHS of Eq. 



-8- 



may cancel each other out) to ensure that the zero is not associated with a curvature 
singularity. 

In the case where 7 7^ 0, we can also write, 

^(5) ^ e2a(r,)^(4) ^ M (20 {a'iv))^ - 8a" iv)) (4.6) 
= -|7l(ymr/))+(l>'(r/))'). (4.7) 

5. Partially analytic solution with dS4 slices and no curvature singularities 

Having demonstrated above that the dS4 case must produce metric zeroes (or a divergent 
metric) in the bulk, to proceed with our analysis of that case we now no longer take the 
warp factor to be the exponential of a real-valued function. Instead, we work directly with 



the warp factor uj introduced previously, employing Eqs. |2.27| - p.28| as the equations of 
motion. 

Let us denote the roots of w(r/) by r/o, i.e. uj{rio) = 0. Then in the dS4 case where 7 is 
positive, it follows immediately from Eq. |2.27| that 

c^'(r?o)=±l, (5.1) 

for the case of interest where lo" and remain finite at the metric zeroes. Differentiating 
both sides of the equation of motion up to fourth order and evaluating at y = yo, we obtain 
the relations: 

c^"(7?o) = 0, (5.2) 
l>'(7?o) = 0, (5.3) 
io^^Hvo) = 0. (5.4) 



(See also for an alternate approach yielding Eq. |5.3| ). 

If ^*(r/) is not 1:1, then it cannot be inverted to yield an expression for y as a function 
of and thus V may fail to be a single- valued function of It is still possible that V might 
be a single-valued function of $ in any case, but we can ensure this is true by requiring <I> 
to be nondecreasing. With this additional constraint, any point where $'(r/) = must also 
be a turning point in the derivative, i.e. *&"(??) must also be zero. 

This assumption has consequences for the form of the required potential. We can 
rewrite the Klein-Gordon equation (Eq. ^.24 ) in terms of w(r/). 



3^ = i"(„)+45:^j,.(,). (5,5) 

If $(r/) is nondecreasing and continuous, then <I>'(ry) > OVry, and as zeroes in the metric 
correspond to points where uj'{rj)/uj{i]) changes sign, it follows that zeroes in the metric are 
either zeroes in the derivative of the potential with respect to ^> (if dF/d^> is continuous 
at these points), or cusps where dy/d$ changes sign. This latter case was proposed as a 
possibility by Davidson and Mannheim [p6|. 



-9- 



Imposing the condition that $"(t/o) = and differentiating the equation of motion up 
to sixth order yields the further relations: 

[j'Hrjo))' = J'\vo)co^'Hvo), (5.6) 

J''Hm) = o. (5.7) 

Note that the conditions on the metric are a result of our somewhat artificial approach 
in starting from the metric and attempting to reconstruct a physically reasonable scalar 
field kink and potential. It is not clear to what degree they indicate constraints on the class 
of potentials that can give rise to a scalar kink supporting a dS4 domain wall. However, 
we have also derived constraints that apply directly to the potential, implying that it must 
have turning points or points of inflection at the zeroes of the metric. This constitutes a 
limitation on the physical potentials that can give rise to a smooth singularity-free scalar 
field and metric with dS4 cosmology on the brane. In flat space we would expect V{^) 
to possess minima at the asymptotic values of the scalar field ^(rj), that is at the limits 
lim^_^-l-oo ^(??)i but, as noted above, in the present case the potential must have extrema 
at the zeroes of the metric. The values of $(??) at the zeroes of the metric will generally 
not coincide with the asymptotes, so a simple double- well potential cannot give rise to the 
desired metric and scalar field configuration. 

A straightforward but inelegant method for writing down a warp factor satisfying these 
conditions is to start with a simple cosh function, motivated by the warp factor in the case 
of an infinitely thin brane with a time-dependent metric (P, |^, |3^). Adding additional 
linearly independent even functions (sech, sech^ and sech^ functions are suitable) which 
decrease rapidly far from the brane and are zero at the roots of the warp factor, we retain 
the essential features of the delta- function brane solution, but the coefficients of these terms 
can be adjusted to fine-tune the derivatives of the metric as required. 



Once the metric has been written down, Eq. 2.27 immediately yields the derivative of 
the dimensionless scalar field ^'{rj), and integrating gives ^{r]) itself. Equation 2.28| then 
yields the dimensionless potential F as a function of rj, and inverting the function <I*(r/) 
yields an expression for V as a function of 

However, at this stage the integration to obtain ^{ij) (and consequently the inversion 
of this function to obtain V{^)) has not been performed analytically, due to the inelegance 
of our method for writing down a metric satisfying all the requirements of the problem. 

Specifically, the metric ansatz, 

uj{r]) = A — Rcosh{rri) + B ^sech(r77) — — ^ + C ^sech^(r77) — (^~^ ^ + 

+D ( sech^(rr?) - ( ^ ) \+E[ sech'^(r7?) - ( ^ ) ) , (5.8) 





satisfies Eqs. pl - [5.2| and 5.4-5.6, with A and R as adjustable parameters and B, C and D 
given by, 

B = {A^ (16^1° - 2mR^A^ + UMR'^A^ - 3327i2^A^ + mUR^A^ - 1400i?^°) + 
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+ER^ (-8192^^2 + 73728A^°i?2 - 279040^^ i?^ + 567296 A'^i?'^- 
-649152^^ i?8 + +393792^2 - 98448i?^2^} / [A^R (16^^^ _ iqsA^^R^+ 
+766^^ i?^ - 2019A^i?^ + m^A^R^ - 2432^2 + 768i?^2)] , (5.9) 
C = ]^[A^ [29, A^ - A4:m^A^ + Uim'^A^ - lUm^A^ + 784i?'^) + ER'^ (5600y4^2_ 

-53648y4^°i?2 + 217S2SA^R^ - 479506^^ i?^ + 595020A^i?'^ - 2>m42AA^R^^+ 
+104160i?^2^} / [A^ (16^12 _ i68yli0i?2 + 7mA^R'^ - 2019^*^ i?^ + 3070A^i?^- 
-2432^2^^0 ^ 768^i2)j ^ (510) 

D = 1 {A^ (2A^ + im^A^ - mR'^A* + 128i?M2 - 64/?^) + Ei?^ (-896^^^^ 

+9056>l^°i?^ - 39864A^i?^ + 98164y4^i?^ - 137328y4'^i?^ + imAlQA^R^^- 

-29568i?^2)} / [R^A^ (16^^^ - l&SA^^ R^ + imA^R^ - 2019^^ i?^ + m^A^R^- 
-2432^2^^0 ^ 768^12^] _ (511) 

The parameter E is obtained by solving the quadratic equation, 

= X {MdSR^'^A^ - 11424i?^^) + 

+E X (-1225i?^A^^ + bllQBR^A^ - 2MSCR^A^ - Sll2CR^^A^ + 5488^^ i?*^ - 
-5040A'^i?^ - 1225Si?^^^ + 12832L»E^°A'^ - imdDR^A^ - A^BR^A^ - 
-lUmDR'^'^A^ + m^QCR^A^) + 

+lQmD'^R^A^ + SB'^R^A^'^ - SBRA^^C + nA^^CR^ + lOOOA^^Di?'^ - 
-2QA^'^BR^ - 112{)A^DR^ - 3C^^^ - lb{)DR^A^^ - UQA^^CR^ + 
+7QBR^A^^C - 2080Z?^i?^°A^ - UAC'^R^A^ + SA^^BR - ^mBR'A'D - 
~112BR^A^C + 912BR^A^D - IbOBR^A^^D + 1568CR'^A^D + 96C'^R'^A^'^ - 
~1760CR^A^D - 192CR^A^^D - UB^R^A^^. (5.12) 

The parameter r is then determined by solving the equation, 

, ,,2 r^A^ -R^)(A'^ + BRA^ + 2CR^A^ + ADR'^A'^ + 6ER^f , , 
1 = (w (r/o)j = . (5.13) 

For example, setting ^ = 5, i? = 1, we obtain, 

B K. 24.0, C « 2.53, D « 8.17, E ^ -1.66, r « 0.101. (5.14) 

Figures ^ and |2| illustrate the bulk profile of the metric and dimensionless scalar field in 
this case, while Fig. ^ depicts the required potential. 

6. "Locally localised gravity" solutions with AdS4 and dS4 brane cosmol- 
ogy 

In the case where 7 < and the brane cosmology is anti-deSitter, the metric zeroes which 
complicate the dS4 case do not occur and it is possible to write down a fully analytic solution 
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system. The simplest such solution (which may also be obtained by a superpotential 
approach starting with a trigonometric periodic potential) has the warp factor exponent 
(T(r/) = A — lncosh(rr/). However, we shall first consider a slightly more complicated trial 
warp factor, and recover this simple solution as a special case. 

Karch and Randall ( [^] ) discussed warp factors giving rise to localised four-dimensional 
gravity in the case of an infinitely thin brane. In the case of a brane with AdS4 cosmology, 
the warp factor grew exponentially far from the brane, like the simple solution described 
above, or like the warped metric presented in [31| generated by a "ghost" scalar field. How- 
ever, close to the brane the metric behaved qualitatively like the decreasing warped metric 
associated with a Minkowski brane (Sec. ^, [Q]), and it was this local behaviour that was 
responsible for confining gravity to the brane. It seems reasonable, then, that a similar but 
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Figure 3: Potential V{^) generating dS4 brane cosmology; A = 5, R = I, see Eqs. 5.8-5.14 



smooth metric (with a decreasing warp factor close to the brane that then turns around 
and increases exponentially) might be required to localise four-dimensional gravity on a 
scalar field domain wall in the AdS4 case. An example of such a warp factor exponent is 
plotted in Fig. and graviton confinement in coupled metric-scalar systems of this form 
was studied by Kobayashi et al p^]. 

It is possible to "smooth out" warp factors initially derived for a thin-brane system, 
essentially by replacing \rj\ with ln(2 cosh rj). Using this idea as motivation, we can construct 
exact analytic solutions for the scalar field and metric with AdS4 and dS4 slices: the 
supporting scalar field is kink-like in the AdS4 case but contains curvature singularities in 
the dS. discussed above. 

We consider a trial warp factor of the form, 



uj{r]) = Acosh.{rr]) + Bsech.{rr]). 



(6.1) 



IfO < A/B < 1, then this metric has off-brane turning points similar to those in the 
solutions of |3^, (see Fig. The limit where B = corresponds to the simple solution 



mentioned at the start of this section. In any case, Eq. 2.27 becomes: 

2 



1 



1 



(Acosh(r?7) -|- Bsec\i{rri)) 

X 



1^ {-A{A + AB) + AABsedi^irr]) + B'^sech'^{rr])) 



(6.2) 
(6.3) 



To ensure that the RHS of the equation is always positive and to facilitate an analytic 
solution, we impose the relation 



^=r^ UA^ + A(A + AB)) 



(6.4) 
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Eq. 6.3 then yields, 



(*'(>)) 



V A + Bsech^{rr]) J 



(6.5) 



and taking the square root of both sides and integrating yields the scalar field profile, 
arctan sinh(rr/) + \ — — arctan [ \ — — sinh(r?7) 



A + B 



A + B 



(6.6) 



The dimensionless potential V can easily be written as a function of r/, 
/ \ 2 



[A + Bsedi^{rr])) 



2A^ - 2A{2,A + 2B)sech2 (rr/) 



2B{A + B)sech'^{rr]) + -B'^sech^{rr]) 



(6.7) 



however ^*(r/) cannot easily be inverted analytically, so V can generally not be written 
analytically in terms of standard functions of $ (although $ is 1 : 1 and therefore always 
invertible, so V can always be expressed numerically as a well-defined function of $). 

AdS4 (dS4) brane cosmology corresponds to the case where 7 < (7 > 0), and by Eq. 



6.4| , this is equivalent to requiring, 

+ iAB > 0, AdS case, 
+ iAB < 0, dS case. 



(6.8) 
(6.9) 



In the AdS case, Eq. |6^ can obviously be easily satisfied provided A,B>0. In this 
case ^y A/ {A + B) is real, and Eq. |6.6| yields a kink-like profile (Fig. |5|). This profile can 
be numerically inverted to give a plot of the generating potential V{^) as a function of 
the scalar field $ (Fig. P). As noted previously, the case < A < B yields a metric 
qualitatively similar to that of [^], which is expected to localise gravity |22]. 

In the special case where B = 0, Eq. |6.6| becomes, 

(6.10) 
by noting that 

(6.11) 



^(rj) = ±2 arctan sinh(r?7). 
Then the potential V{^{y)) can be obtained analytically from Eq. 

1 



sech [ri]) 



cos 



1 + sinh^(r7/) 

and so in this special case, the dimensionless potential takes the form. 



^(77) 



-2r' 



l + 3cos2 (I 



(6.12) 



In fact, it is possible to obtain a slightly more general solution in this case: if we do 
not require Eq. to hold, then the scalar field configuration for i? = 0, 7 < becomes, 

±Vl - A^r'^ 



$(ry) 



Ar 



■ arctan sinh(rry), 



(6.13) 
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Figure 4: Warp factor exponent 17(77) for a domain wall with AdS4 brane cosmology; A = 0.2, B = 
0.8, see Eq. O 



2- 
1- 




-10 -8 -6 -4 -2 "/ 

/-I- 

y -2- 


2 4 6 8 10 



Figure 5: Dimensionless scalar field profile <&(??) for a domain wall with AdS4 brane cosmology; 
A = 0.2, B = 0.8, see Eqs. 



and Eq. 2.22 yields, 



y(<l) 



-2r" 



cos 



AR 



Vl - A2r2 



(6.14) 



which reduces to Eq. |6.12 in the case where bA^r"^ = 1. Note that in these cases the 



asymptotes of the scalar field do not coincide with local minima of the potential. 



Now in the dS case (7 > 0), it follows from Eq. [6]^ that 4A{A + B) < 0, and 
consequently A/{A + B) < 0. A sample profile of the warp factor exponent is shown in 
Fig. 1^: as expected from the results of Sec. it diverges to cxd a finite distance from the 
brane. (The warp factor exponent is obtained by taking — lna;(r/); we plot the exponent 
rather than the warp factor itself to facilitate comparison with previous studies.) 
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Figure 6: y(3>) vs $ for a domain wall with AdS4 brane cosmology; A — 0.2, B — 0.8, see Eqs. 



6.1-6.6 
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Figure 7: Warp factor exponent a(r]) for a brane with dS4 cosmology lying between two naked 



singularities; A — ~1, B — 2, see Eq. 6.1 



Writing ia = A/(A + B), for a real, the scalar field profile of Eq. S.6 becomes, 

$(7/) = lb [arctan sinh(r?7) — a arctanh (a sinh(rr/))] . (6.15) 

Let us assume a + sign for the purpose of this analysis: there are no significant differences 
between the kink and antikink solutions. A sample profile is given in Fig. ^. 

As is clear from Fig. |8|, the second term diverges to ±00 where asinh(r77) = ±1: 
it is trivial to check that the points of divergence correspond precisely with the zeroes 
of the metric. Close to these singular points, which we shall denote ±r/o, can be 

approximated by, 

l>(-?7o + Ar/) = - arctan ^-^ ~ f Qra\/l + a^Ar/^ + 0(Ar/), (6.16) 
4>(r/o - Ary) = arctan ^- j + ^ In Qra\/l + a^Ar?^ + 0{Ar]). (6.17) 
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Figure 8: Dimensionless scalar field profile ^(r/) supporting a brane witli dS4 cosmology; A 



-l,B^ 2, see Eqs. 6.1, 6.15 



Inverting these expressions allows the potential V to be written analytically as a function 
of close to the singularities. We find that the leading order term in V{^) is given by, 



arctanf — ) ^ — ll*! 



(6.18) 



so in this case t^(^) is unbounded below, and the scalar field diverges a finite distance from 
the brane. The potential can also be plotted numerically against the scalar field over the 
range between the singularities (Fig. lol), but there are no other features of real interest. 
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Figure 9: V^(^) vs <f> for a potential supporting a brane with dS4 cosmology; A = —1, B = 2, see 

Eq. [elsl . 



The 5D curvature scalar is given by Eq. 4.7, and close to the singularities the leading 



order term is of the form, 

^^'^(^)~l7l«'-r^. (6.19) 
(Ar/) 

Thus in this case the zeroes in the metric represent curvature singularities, and indicate 
divergences in the scalar field energy-momentum tensor, rather than merely being horizons. 

7. Fermion trapping by warped metrics 

The trapping of fermions on a scalar field domain wall was first described in the 1980s (||l|, 



32]). More recently, fermion trapping in a scalar field domain wall with a five-dimensional 
warped metric has been discussed by a number of authors (||l^, |l8|, 3C, ^]), but to our 
knowledge only in the case of a time-independent metric. 

To write down the Dirac Lagrangian in curved spacetime, we employ the vielbeins 
and inverse vielbeins V^^ defined by (|3^, |39[|): 

QMNix) = Fi(x)V^(x)r/AB. (7.1) 

Here ijab is the Minkowski space metric, and A,B,C,... indicate coordinates in (locally 
defined) 5D Minkowski space, with M,N,... indicating coordinates in the curved space 
described by the metric guN- 

The Dirac Lagrangian can then be written out as, 

= ^r^VA'if - gF^F{^)^, (7.2) 

where the F's are the Minkowski space Dirac matrices, F is some odd function of ^ de- 
scribing the coupling between the Dirac field and the scalar field, and T>a is the covariant 
derivative with spin connection defined by, 



^A^Vr^^ + ^M{x)], (7.3) 
1 



Em(x) = ^a^^<ycJV;M, (7.4) 



VcN;M = ~ ^NM^CR-, (7.5) 



J- NM — 2^ 



dgMN _|_ dgNS _ dguN 
dx^ dx^^ dx^ 



(7.6) 



Here describes how the field transforms under infinitesimal Lorentz transformations, 
i.e. the spin of the field. For a spin-1/2 field, cj^^ = (l/4)[r^, T^]. 

We may choose the Dirac matrices in five-dimensional Minkowski space, F^, to be: 

r" = 7", = (7.7) 

Here the 7 matrices are simply the 4-dimensional Dirac matrices, which obey the four- 
dimensional Clifford algebra {7°^, 7^} = 2r/"^ with our chosen metric signature ( — h ++). 

It is easily verified that these matrices obey the five-dimensional Clifford algebra, 
|pA ^ 2r]^^, and it follows that for a Dirac field, ii A ^ B ^ a^^ = (l/2)r^r^ (for 
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A = B, of course, a"^^ = 0). Note also that raising or lowering the indices on the F's is 
trivial: = -Tq, = Ta for Aj^O. 

As previously, we consider a 5D metric of the form given in Eq. 2.1. For now, we shall 
employ the o"(y) parameterisation, giving the key equations rewritten in terms of at 
the end of this section. Let us choose a set of 4D vielbeins satisfying, 

(7.8) 

Then the 5D vielbeins and inverse vielbeins may be chosen as, 

ij' ~\ 1 

With this choice of vielbeins, the five-dimensional spin connection T,m becomes: 

+ icT'(2/)e-(2')r4r°t>„^ 





„(4) _ a 

yfii/ ~ ^/^ '^ly 'lap- 



-'M 



(7.9) 



(7.10) 



where S„ is the spin connection for fermion fields in the four-dimensional spacetime de- 



scribed by the metric ^/j^ 

Consequently, the covariant derivative with spin connection becomes 



d_ 



and the sum T^Da can be written. 



T^Va 



7"< 



+ 75(^-2cT'(y)). 



(7.11) 



(7.12) 



The five-dimensional Dirac equation for a fermion field coupled to gravity and the 
background scalar field ^ is simply, 

[F^Va - 9fF{^)) ^ = 0. (7.13) 

It is well known ( |lj, 0, ^] ) that the confinement of fermion modes to the brane 
depends on their chirality. Let us therefore write the Dirac spinor ^ in the form 

^{x)=UL{y)Mt,oc') +UR{y)i;R{t,x'), (7.14) 

where ipi and ipR are the left-handed and right-handed components of a four-dimensional 
Dirac field, and hence j^^l = —"^L, TsV'-R = V'-R- Inserting this ansatz and Eq. 7.12| into 
the Dirac equation yields, 

= {T^Va - 9fF{<^)) ^ 



7 






d 




) 










d 





ll^Lit,x' 



lpRit,x') + 



+^pR{t,x') ^■^-2a'iy)jUR{y)- 
-gFF{^)UL{y)Mt,x') - gFF{<^)UR{y)ijR{t,x'). 



(7.15) 
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Now suppose we require that the four-dimensional spinors 'i/'Lj '4'R satisfy the Dirac 
equation for fermions in the four-dimensional spacetime described by the metric g^/^J , 

7"< {d^ + S^) Vl = mVi?, (7.16) 
7X(5A. + S^)V'i? = "iVl- (7.17) 
Then the five-dimensional Dirac equation becomes, 



= tPRit,x') 



me'^^y^URiy) " " 2^7' (y)) U^y) - gFF{'^)UL{y) 



(7.18) 



and equating the coefficients of ipR and V'L separately to zero, we obtain a pair of coupled 
first-order differential equations in the bulk coordinate, 

me'^^y^ULiy) + - 2t7'(y)) i^dy) - gFF{^)UR{y) = , (7.19) 
^"^^URiy) - - 2a' (y)^ UL{y) - gpFmUUv) = . (7.20) 



me" 



Note that these equations are identical to those for a static domain wall (see for example 
1 14]): the generalisation to any 4D metric g^^y only modifies the warp factor cj(y) and the 
scalar field 

Let us now define the dimensionless rescaled mass and coupling constant for the case 
7 7^ 0, in the case where F{^) = 

/ 3 'igp 



Then Eq. 7.1£- |7.2d| can be rewritten in terms of dimensionless quantities, 

^e'^^'^^ULiv) + - 2f^'('?)) ^Riv) - hpHvpRiv) = 0, (7.22) 

//e^(^)Z^ij(7?) - - 2a' (r/)^ Uiiri) - hpMvWv) = 0. (7.23) 

As previously, in configurations containing bulk singularities we may wish to rewrite 
the metric in terms of the warp factor to: the only reason for casting our equations in terms 
of the warp factor exponent is for easy comparison with the existing literature. In this case 
the mode equations can be derived as, 

^^L{y) +(^ + 2441 ^Riy) - gpF^URiy) = , (7.24) 



uj{y) \dy uj{y) 

my) -(^ + 2^4 ) my) - gpF^ULiy) = o , (7.25) 



u){y) \dy uj{y) 

or in terms of dimensionless quantities for the F[^) = ^ case, 

^ ^^l{v) + (^ + 2^1 mv) - hF^v)mv) = 0, (7.26) 



^(ry) \dri ^^i^l) 

^ Mv) -(^ + 2^ ) mv) - hFHv)mv) = 0. (7.27) 



uj{ri) \dr] w(ry) 
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8. Confinement of the fermion zero mode on AdS4 and dS4 branes 



In the case where m = 0, the equations for Ul and Ur simplify by decoupUng, 



^R{y) 
dy 



Uiiy) {2a'{y)-gFF{^{y))), 
UR{y) {2a'{y)+gFF{<^{y))), 



allowing us to analytically study the confinement of the fermion zero mode. We are pri- 
marily interested in the behaviour of the fermion field within the region close to the brane 
and, in the case of dS brane cosmology, bounded by the singularities. We shall thus use the 
a parameterisation, and demonstrate that the presence of even the coordinate singularities 
discussed in Sec. ^ leads to the failure of the usual confinement mechanism. 

The above first-order linear differential equations can easily be solved for the bulk 
coefficient functions, 



(8.3) 



If we substitute these expressions back into the kinetic-energy term of the Dirac La- 
grangian, we obtain, 



7°< 



lPR{t,x') + 



+gFFi<^)ALARe'''^y'> {ijR{t, x')Mt, x') + Mt, x')i^R{t, x')) . (8.4) 
The action for the fermion field is simply. 



d^xy^-g{x)Cii,{x). 



(8.5) 



U g = Det{gl^J), then g = e'^'^'^y^g, and = e"^'^^^)^/^. Noting that the term in 

-F($(y)) is odd for a kink solution, we need only consider the first two terms in the action, 
because the third integrates to zero. These terms factorise into the usual 4D action over 
the brane coordinates. 



d x\f^-g 



V'L/R(i,^*)(7X 



d 



multiplied by integrals over the bulk coordinates given by 



dy e' 



r(y) -2gf./F($(j/))dy 



42 



dy e 



(8.6) 



(8.7) 



for left-handed and right-handed fermion fields respectively. 
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For the case where = $ and 7 7^ 0, the normahsation integrals can be written in 

terms of the dimensionless quantities, 





It is now clear that in general, zeroes in the metric (points where cr{r]) diverges to 
00) correspond to singularities in the integrand, and thus the kinetic energy part of the 
action may be expected to be non-normalisable in cases where the metric contains zeroes. 



In particular, this will be the case for the new solution presented in Eqs. |5.8| - p.l4| . This 
obviously represents a challenge for that kind of model: although it has the nice properties 
that the bulk inverse-metric singularity is a coordinate not a curvature singularity and the 
scalar field configuration is smooth, it has the drawback that the usual kink-based fermion 
zero-mode localisation mechanism does not generalise to that dS4 case. Dealing with this 
challenge is beyond the scope of this paper, but two logical approaches immediately suggest 
themselves: One could look to introduce new physics, beyond the scalar kink, to localise 
fermions. Alternatively, we know that our universe is only approximately de Sitter, so it 
would be interesting to explore fermion localisation via the kink mechanism in an effective 
domain wall FRW cosmology that has the transitions from radiation to matter to vacuum 
energy domination of regular cosmology. 

It seems possible that if the kink <I> also diverges to ±00 at the metric zeroes, then 
one of the chiral fermion zero modes might be normalisable over the bulk. However, in 
the dS4 solution we have described with curvature singularities at the metric zeroes, this 
cancellation does not occur. Eqs. |6.16[]6.17| show that close to the singularities, the leading 
order term of ^{ij) is of the form ln{qAr]), so the integral function J $(r/) drj is bounded 
close to the singularities, while e"^^^^ diverges as 1/Ar/. 

If the warp factor exponent cr{y) is smooth and nonsingular, the fermion zero modes 
may or may not be normalisable. Suppose ~ c\y\ as y ^ ±00, for some c < 0. Then 
one of the two chiral fermion zero modes will always be confined: the chirality of this 
mode will depend on the scalar field profile. The other mode may or may not be confined, 
depending on the asymptotic behaviour of J dy F{^{y)) relative to <T(y). This is the case 
when the metric diverges to infinity asymptotically, and is relevant for the cases considered 
in Sec. with AdS4 cosmology on the brane. This behaviour has also been studied by 
Koley and Kar in the context of a "ghost" scalar field with negative energy in the bulk in 



1 31], and by Bajc and Gabadadze for fermions localised on a non-fine-tuned RS brane by 



a scalar field in [33|. 

If conversely the warp factor exponent behaves as (T{y) ~ c\y\ for y — > ±00, with c > 0, 
then one of the two modes will certainly be non-normalisable, while the other may be 
normalisable depending on the behaviour of the integral J dy F{^{y)). This is the case for 
fine-tuned static brane solutions ([0, |3^)j where the metric elements approach zero 
for large \y\. This type of warp factor is employed to localise gravity on the brane in the 
Randall-Sundrum approach (0). 

In particular, if F{^{y)) has a kink or anti-kink profile, and limy^oo F{^{y)) = 
— limy^^oo F{^{y)) = A, then far from the brane, J dyF{^{y)) ~ X\y\ (up to a con- 
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stant of integration). Consequently, if cr{y) ~ c\y\ for \y\ large, then the normalisation 
integrands are of the form, 

exp{{c + 2gFX)\y\), exp ((c - 25fA) |y|) , (8.9) 

for right- and left-handed fields respectively. 

For the "locally localised gravity" AdS4 solution outlined in Sec. |^, these conditions 
hold for the simple coupling -F(<&) = ^, with the parameters, 

\/3 TT , . . . 

A = — -(1 + a), c = -r. (8.10) 
Thus the action is normalisable for right-handed fermion fields provided, 

V3vr(l + a) 

and for left-handed fermion fields if 

V3vr(l + Qj 

but the latter relation is always true for positive gp. 

In terms of the dimensionless quantities, we replace y with ry and gp with hp m the 
discussion above, and the parameters become c = — r, A = 7r/2(l + a). The normalisability 
conditions become, 

hp < —TT- T (8-13) 



for right-handed fermions, and 



hp > -— -, (8.14) 

vr (1 + aj 

for left-handed fermions. 

However, normalisability of the chiral fermion zero-modes does not necessarily imply 
that those modes are localised on the brane. For example, in some parameter regimes, the 
normalisation integrand has two peaks at the points where the warp factor turns around, 
and a local minimum (rather than a local maximum) on the brane, before falling off rapidly 
far from the brane. If the local minimum is sufficiently shallow and the peaks are sufficiently 
close to the brane, this behaviour may still constitute localisation to the brane. 

The Lagrangian density always has a local extremum at the brane (i.e. its derivative 
with respect to the bulk coordinate is zero there): its second derivative with respect to the 
bulk coordinate determines whether the brane lies at a local minimum or maximum. The 
second derivative of the integrand of Eq. |8.8| has the same sign as 

a"{rj)±2hp^'{rj), (8.15) 

where as previously, the + sign applies to right-handed fermion fields and the — sign 
corresponds to left-handed fields. 

For the AdS4 case outlined in Sec. ^, 

(t"(0) ± 2hF^'{0) = (1 - 2a2) ± 2hpr (l + a^) , (8.16) 
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so the brane lies at a local maximum of the Lagrangian density for, 



±hF < 



r {2a^ - 1) 
2(a2 + i) 



.17) 



Note that in the case where < A < B, a < 1/2 and therefore the RHS of this expression is 
always negative. Consequently, this condition cannot hold for right-handed fermion modes: 
this result is analogous to the usual behaviour of fermion zero modes for a Minkowski 
brane, where right-handed modes are always unconfined (fl^, |3^). Fig. 10 demonstrates 



a sample profile of the integrand of Eq. |8.8| for the case where the right-handed fermion 
modes are normalisable, but not confined to the brane. Left-handed fermion modes are 



always normalisable, and may (Fig. [Tl| ) or may not (Fig. 12) exhibit a peak at the brane 
itself. 




Figure 10: Bulk dependence of the kinetic energy term, Eq. 8.8, for right-handed fermions on a 



domain wall with AdS4 brane cosmology; A = 0.2, B = 0.8, hp = 0.1; see Eqs. 6.1-6.7 and Eq 



9. Conclusion 

We have presented solutions to the Einstein equations describing a scalar field coupled to 
five-dimensional gravity, with a warped five-dimensional metric and dS4 and AdS4 brane 
cosmology. In the dS4 case, the metric necessarily contains zeroes which generally (but 
not inevitably) correspond to curvature singularities. It is possible to obtain a warped 
metric with dS4 slices from a smooth potential and scalar field, and in this case the metric 
zeroes are simply horizons. However, there are stringent conditions on the metric in this 
case which make it difficult to write down an analytic solution. Moreover, the presence of 
zeroes in the five-dimensional warped metric, whether representing horizons or curvature 
singularities, generally leads to divergences in the normalisation integrals for both the left- 
and right-handed zero modes of the Dirac field. This issue does not arise in the case of an 
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Figure 11: Bulk dependence of the kinetic energy term, Eq. B 
domain wall with AdS4 brane cosmology; A = 0.2, B = 0.8, hp 



for left-handed fermions on a 



0.3; see Eqs. 6.1-6.7 and Eq 



7.21 




for left-handed fermions on a 



Figure 12: Bulk dependence of the kinetic energy term, Eq. 
domain wall with AdS4 brane cosmology; A = 0.2, B = 0.8, hp = 0.1; see Eqs. 6.1-6.7 and Eq 
7^. 



infinitely thin brane, where the matter fields are confined a priori to a 3-|-l-dimensional 
slice of the higher-dimensional space and do not extend into the bulk at all. 

In the AdS4 case, we have derived new analytic self-consistent solutions for the metric 
warp factor and scalar field kink, motivated by the principle that localisation of gravity 
should depend only on the behaviour of the metric close to the brane (||3^). Although 
for non-trivial brane cosmology the scalar field does not generally seem to asymptote to 
minima of the potential, these configurations (and their dS4 counterparts) are expected to 
be both classically stable and to confine the four-dimensional graviton p2| . The ensuing 
warped metric admits both left- and right-handed normalisable fermionic zero modes (as 



in ^]), although not all the normalisable modes are localised on the brane. At present 
we have only investigated massless chiral fermions, but it might also be interesting to 
investigate the spectrum of massive fermionic modes in these backgrounds, using Eqs. 



7.1C-7.2C. 



Our present analysis only applies to a particular class of warped metrics, corresponding 
to four-dimensional metrics of constant curvature. Future studies dealing with the effects 
of matter on the brane will need to employ a more general metric ansatz, allowing metric 
elements that are non-separable functions of the brane and bulk coordinates. 
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